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The enhancement of tunneling probability in the nearly integrable system is closely examined, 
focusing on tunneling splittings plotted as a function of the inverse of the Planck’s constant. On 
the basis of the analysis using the absorber which efficiently suppresses the coupling creating spikes 
in the plot, we found that the splitting curve should be viewed as the staircase-shaped skeleton 
accompanied by spikes. We further introduce renormalized integrable Hamiltonians, and explore 
the origin of such a staircase structure by investigating the nature of eigenfunctions closely. It 
is found that the origin of the staircase structure could trace back to the anomalous structure in 
tunneling tail which manifests itself in the representation using renormalized action bases. This also 
explains the reason why the staircase does not appear in the completely integrable system. 

PACS numbers: 05.45.Mt,05.45.-a,03.65.Xp,03.65.Sq 


I. INTRODUCTION 

The tunneling effect is peculiar to quantum mechanics 
and no counterparts exist in classical mechanics. Quan¬ 
tum tunneling plays a role and actually manifests in vari¬ 
ous situations ranging from atomic and molecular physics 
to various phenomena in condensed phases. 

In most cases, incorporating the tunneling effect into 
each case is made by using the system with a single degree 
of freedom. This is justified and certainly provides a good 
description if the tunneling penetration proceeds in only 
one direction, but this is not the case when the system 
has multi degrees of freedom. 

The most important qualitative difference between 
one- and multi-dimensional systems would be that classi¬ 
cal particles are confined not only by the energy barrier, 
but also by the dynamical harrier. The latter is formed 
when additional constants of motion, either globally or 
locally, exist besides the energy. What is more crucial 
is the fact that generic multi-dimensional systems are 
no more completely integrable and chaos appears in the 
underlying classical dynamics, so one must take into ac¬ 
count new aspects of quantum tunneling absent in com¬ 
pletely integrable systems [l, ,2|. 

The role of classical chaos in quantum tunneling has 
first been discussed in the observation of the wave packet 
dynamics Q, and then clearly recognized in the behavior 
of the tunneling splitting of eigenenergies |J-[6l]. To un¬ 
derstand why chaos could play a role in the tunneling pro¬ 
cess, it suffices to suppose the states forming a doublet, 
which is a complete analog of the doublet appearing in 
the system with one-dimensional symmetric double well 
potential. It is important to note, however, that the dou¬ 
blet in multi-dimensional systems are supported by sym¬ 
metric regular tori in phase space and chaos exists in be¬ 
tween. As one varies an external parameter of the system, 
it can happen that states forming the doublet and a state 
supported by the chaotic region come close to each other 
in the energy space and form avoided crossing. Within 
the interaction regime, the energy splitting between the 


doublet becomes large through couplings with the chaotic 
state, meaning that the tunneling amplitude between one 
torus to the other is enhanced. Chaos-assisted tunneling 
(CAT) occurs in this way [ 5 , 5[. A similar enhancement 
is known to take place if the doublet is bridged by nonlin¬ 
ear resonances. Nonlinear resonances are also important 
ingredients in multi-dimensional systems, and the latter 
mechanism is called resonance-assisted tunnelinq (RAT) 

Effl. 

In order to go beyond qualitative explanations and 
to obtain more direct evidence for the connection be¬ 
tween chaos (and/or nonlinear resonances) semiclassical 
(WKB) analyses are desired, especially based on the com¬ 
plex classical dynamics since quantum tunneling is a clas¬ 
sically forbidden process. 

For the system with one degree of freedom, there in¬ 
deed exists a standard approach that has been estab¬ 
lished already [13]. The instanton is the name of a com¬ 
plex orbit which conveys the tunneling amplitude run¬ 
ning along the imaginary time axis, and the formula rep¬ 
resenting the tunneling splitting in the symmetric double 
well mode is expressed as 

A E ~ ahe ~ s/h , (1) 

where and classical ingredients a and S can be deduced 
in the instanton calculation (l3j . 

On the other hand, in multi-dimensional cases, full 
semiclassical analyses using the complex classical dynam¬ 
ics could so far be applied only to the time domain [1411151 ] 
and have not been even formulated except for completely 
integrable situations in the energy domain [lG - ilOf . The 
enhancement of tunneling could therefore be well ac¬ 
counted for in terms of the complex dynamics [20j . but 
fully convincing semiclassical understanding for the en¬ 
ergy domain is still lacking. 

The aim of this article is to explore the origin of the 
enhancement of tunneling probability observed in the en¬ 
ergy domain. The tunneling probability in the energy do¬ 
main is often measured by tracking the energy splitting 
or the properly defined tunneling rate as a function of 
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1/h. The enhancement is typically observed as plateaus 
accompanied by spikes due to energy resonance [f| [gj]. 
Characteristics observed there are understood within a 
framework of RAT, and in particular, plateaus could be 
interpreted as a kind of phenomena that might be called 
quantum overlapping resonances; a bunching of spikes, 
each of which is associated with an individual quantum 
resonance, turns out to create plateaus, or a persistent 
long-range interaction of each resonance with other states 

11, ll. 


However, one should recall that the enhancement oc¬ 
curs even when the Planck’s constant is not small enough 
to resolve the chaotic components or nonlinear resonance 
islands. This rather paradoxical behavior has been al¬ 
ready observed in several models [HEil, but its origin 
in such a slightly perturbed regime has never been seri¬ 
ously investigated to the authors’ knowledge. 

We here take a close look at the nature of the en¬ 
hancement of the tunneling probability in such nearly 
integrable regimes by introducing techniques such as ab¬ 
sorbing the individual states involved in avoided crossings 
and decomposing the eigenstates into proper integrable 
bases, as explained in detail below. We especially focus 
not only on the behavior of the energy levels but also 
on the nature of eigenfunctions to elucidate which com¬ 
ponents are mostly responsible for the plateau structure 
formed in the energy splitting vs 1/h plot. 

In conjunction with this, we shall stress the importance 
of observing wavefunctions in the whole range because 
there are various ways to define the “tunneling proba¬ 
bility” , and the nature of tunneling may look different 
depending on how it is defined. Here, for the closed sys¬ 
tem (standard map) the splitting of energy levels will be 
adopted to measure the tunneling probability, whereas 
for the open system, like the Henon map, the probability 
in the asymptotic region is naturally introduced, and the 
decay rate for the absorbed system is sometimes used. 
There would be no legitimate way or one should even say 
that providing a proper definition of the tunneling prob¬ 
ability itself is an issue to be explored in nonintegrable 
systems. Therefore, one should examine more carefully 
the tail of wavefunctions in the whole range before focus¬ 
ing on the amplitude at a certain specific position. 

The present analysis is motivated by a recent work in 
which the mechanism of the instanton-noninstanton (I- 
NI) transition has closely been studied in terms of quan¬ 
tum perturbation theory 2l|, and so spirits and tools 
for analyses are overall common. The term instanton- 
noninstanton (I-NI) is named after the first transition at 
which the deviation from the instanton prediction starts 


The organization of the paper is as follows: In sec¬ 
tion ED we introduce the system studied in this paper, 
and present aspects of the enhancement of the tunneling 
probability by observing the quantum number and 1/h 
dependence of the tunneling probability in our model. 
In section m introducing an absorbing operator, which 
projects out a given set of integrable states, we examine 


0.3 


a, 0.0 


-0.3 

0.3 


a, o.o 

-0.3 

-1.0 -0.5 0.0 0.5 1.0 

q 

FIG. 1: (Color online) Classical phase space for the sym¬ 
metrized standard map / with (a) r = 2/3 and (b) r = 1. 
There are no visible resonance chains for r = 2/3 in the inner 
tours region while chaotic regions around the unstable fixed 
point at (q,p) = (0, 0) and some resonance chains (1: 8, 1: 10, 
1: 12,- • •) become visible for r = 1. They are shown in dark 
blue. The black box put in the upper right corner indicates 
the size of the effective Planck’s constant for h = 1/70. 



which states are responsible for creating spikes typically 
observed in the splitting curve and whether or not the 
staircase structure of the splitting curve appears as a re¬ 
sult of local quantum resonances in the energy space. In 
section El we investigate the nature of eigenstates to 
clarify the mechanism of the enhancement by focusing 
on the local probability amplitude of eigenfunctions and 
the contribution spectrum introduced in [2l| ■ In section 
fVl on the basis of analyses made in section HVl we claim 
that an essential difference of the splitting curve exists 
between integrable and nonintegrable systems. In the fi¬ 
nal section, we summarize and provide outlook especially 
toward our forthcoming papers. 


II. ENHANCEMENT OF TUNNELING 
PROBABILITY 

We consider a quantum system described by the evo¬ 
lution operator in a symmetrized form: 

JJ = e --iF(g)r/2 e -iT(p)T e -|V(9)r/2_ (2) 

The corresponding classical dynamics is given as the sym- 
plectic map / := /V(§) o / T (r) o f v {\) where f v [r) : 
(QiP) 'A (q,P+rV'(q)) and f T {r) : (q,p) ( q+rT'(p),p ) 
are trivial symplectic maps. Here the prime stands for 
the derivative of the function. The classical map / 
corresponds to discretization of the continuous Hamil¬ 
tonian flow for H(q,p) = T(p) + V(q) up to the sec¬ 
ond order of the discrete time step r. Thus, the map 
/ has the integrable (continuous) limit r —> 0, and 
much the same is true on the quantum map ([2|). Here¬ 
after we take the potential function as T(p) — p 2 /2 and 
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V(q) = (fc/47r 2 ) cos(27rg) where k is the strength of the 
perturbation. After rescaling as p K > p/r and fcr 2 = e, 
the classical map / turns out to be the symmetrized stan¬ 
dard map [22|, and the time evolution by the unitary op¬ 
erator U can be interpreted as a single period evolution 
of a 5-functional periodically forcing Hamiltonian with a 
period r. 

In the continuous limit r —> 0, the closed area sur¬ 
rounded by the separatrix is given by S = Vk{2/ 7r) 2 . In 
the following argument, we focus especially on the nearly 
integrable regime and a proper integrable limit will play 
an important role as a reference. In most of situations, 
nonlinearity is controlled by changing the parameter r, 
keeping the parameter fixed as k = ko = 0.7458. 

Figure Q] displays classical phase space for typical 
nearly integrable parameter regions. In the case of 
r = 2/3, classical phase space is predominantly covered 
by regular regions and nonlinear resonance chains are 
not visible in this scale. For r = 1, small chaotic regions 
emerge around an unstable fixed point at (q,p) = (0,0), 
and Poincare-Birkhoff chains induced by nonlinear reso¬ 
nances become visible. Relatively large nonlinear reso¬ 
nances in the inner torus region, which represents libra- 
tional motions in the pendulum Hamiltonian H, are 1:8, 
1 : 10, and 1:12 ones, which are marked in dark blue in 
Fig. [ljb). Below we mainly develop our discussion in the 
case t = 1, but essentially the same argument follows for 
other t cases. 

We numerically solve the eigenvalue problem for the 
unitary operator U 

U\V n ) = u n \V n ), (3) 

under the periodic boundary condition on the region 
(q,p) € (—1,1] x (— 1/2 t, 1/2t]. Let N be the dimen¬ 
sion of the Hilbert space space, then to achieve the peri¬ 
odic boundary condition the relation l/2r x 2 /H= 2irN 
should hold, which yields the relation h = 2/Nr. 

Here u n is expressed as u n = e ~ lE ^ T / h where E n (n = 
0,1,2, •••) are quasi-energies, and |'P„) denote the cor¬ 
responding quasi-eigenstates. Hereafter we focus on the 
doublet states in bounded states supported by the in¬ 
ner torus region, each of which is centered at (q,p) = 
(±1/2,0) and energy splittings between them. 

Quasi-eigenstates | , 3/n) have a symmetry with respect 
to the mirror transformation Tl q : q 1 —> —q, and we 
hereafter denote the doublet states associated with this 
symmetry by ) and the corresponding quasi-energies 
by E±, which form quasi degeneracy. We therefore as¬ 
sign the quantum number n not to an individual quasi¬ 
eigenstate but to each doublet [ 23 }. The states |'P+) 
and |’P~} respectively represent symmetric and anti¬ 
symmetric states. Note that we have additional sym¬ 
metry with respect to the translation T) : q > q ± 1, 
originating from the periodic boundary condition in the 
g-direction. This symmetry does not induce quasi degen¬ 
eracy in energy, but the states belonging to a different 
translational symmetry class do not interact with each 
other even if they have the same mirror symmetry. 



l/h 

FIG. 2: (Color online) (a) The tunneling splitting A E n is 
plotted as a function of the quantum number n, where n is 
labeled not for an individual state but for a doublet (see the 
text). The value of h for each curve is put in the figure. The 
black curves represent the splitting A En in the limit (r —> 0) 
(b) The tunneling splitting AEo for the lowest doublet 
is shown as a function of l/h. The black solid and dotted 
curves correspond to the splitting Ain the limit (r —> 0) 
and the semiclassical prediction ©, respectively. In (b), we 
put the labels (1), (2), ■ • ■, (5) on each characteristic interval: 
(1) first exponential decay (instanton), (2) first plateau, (3) 
second steeply decay, (4) second plateau, and (5) third steeply 
decay regime, respectively. 

In the continuous limit r —» 0, the eigenvalue equation 
for the Hamiltonian H(q,p) = p 1 /2 + V(q) is expressed 
as 

H(q,p)\Jn) = Et {1) \Jt), (4) 

where eigenstates |J^) are in the same symmetry class 
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as the corresponding IT)/). Here the quantum num¬ 
ber n is, as usual, attached in ascending order of eigen¬ 
value En ’^, so | J±) represents the ground state doublet, 
which we will hereafter focus on. 

For the later purpose, we rearrange the quantum num¬ 
ber n for the quasi eigenstates IT)/) such that the overlap 
|(J)/|T)/)| 2 is maximal. This condition, that is one-to- 
one correspondence between |T)/) and | J)/) is fulfilled for 
the values of r used in the present analysis. 

With increase in the value of 1 /h, the tunneling prob¬ 
ability between |T)/), which could be measured by the 
tunneling splitting A E n = E+ — E~, becomes large in 
several orders of magnitude as compared to those pre¬ 
dicted in the continuous limit. The latter is evaluated as 
A En' 1 = En’^ — E n We notice that the overall be¬ 
havior does not depend on the value of the perturbation 
strength r, although the Planck’s cell can resolve chaotic 
regions and nonlinear island chains in the case of r = 1, 
whereas this does not the case at all for r = 2/3. (see 
Fig. [I]and Fig. [3]). 

We illustrate such anomalous enhancement of tunnel¬ 
ing probability in a nearly integrable regime in two ways. 
First, as shown in Fig. 0a), the tunneling splitting A E n 
is plotted as a function of the quantum number n, and 
in Fig. 0b) the splitting AEq as a function of l/h. The 
latter is known to be a standard plot often used in the 
study of RAT (8l4T2| . 

In Fig. 0 a) we notice that, in the relatively large n 
regime, A E n can be fitted by the lines predicted by the 
formula 0, implying that they have completely inte¬ 
grable nature in essence. On the other hand, as n goes 
down from exited states to the ground state, with a fixed 
h , the law described by the formula 0 is violated at cer¬ 
tain critical quantum numbers n c , each of which depends 
on the value of h SE|. 

At such a quantum number n c , the curve for A E n 
changes its slope and forms the plateau. After a certain 
plateau interval, as typically seen in r = 1 and r = 2/3 
for h = 1/80, the slope again becomes large, and then 
forms the second plateau. The emergence of plateaus 
means the enhancement of the tunneling probability as 
compared to the integrable (instanton) prediction. It 
is particularly non-trivial and even paradoxical because 
this enhancement is relatively stronger in the lower dou¬ 
blets than higher excited ones. Note also that the critical 
quantum number n c becomes large with increase in the 
value of r. This sudden departure from integrable tun¬ 
neling has been pointed out in the study of RAT [8], and 
it is called the instanton-noninstanton (I-NI) transition 
in Ref. |2(J [20] , in which the mechanism behind it has 
been investigated in a different perspective. 

The TNI transition is similarly observed in the A E n 
vs l/h plot. As shown in Fig. 0b), the energy split¬ 
ting AEq for the lowest doublet |T[/) exhibits a similar 
behavior. For relatively large values of h , AEq follows 
the instanton prediction Eq. 0, but deviates from it at 
certain values of h, each of which depends on the value of 
r. The staircase-like structure formed with plateau and 


steeply decaying intervals again characterize the overall 
structure. For the purpose of illustration, we call each 
region in the staircase, (1) first exponential decay (in¬ 
stanton) (2) first plateau, (3) second steeply decay, (4) 
second plateau, and (5) third steeply decay regime, re¬ 
spectively (see Fig. 0b)). 

What is prominent in the latter plot than in the for¬ 
mer plot is the appearance of spikes. This is because in 
the former plot Fig. 0a), we could evaluate tunneling 
splitting only at integer values (quantum numbers), so 
may miss spikes even if they exist, whereas we can scan 
AEq at more numerous values of l/h. 

The origin of spikes is a central issue in theory of RAT 
Hi GS GH , in which the effect of nonlinear resonance is 
incorporated by first constructing local integrable pen¬ 
dulum Hamiltonian classically and then applying quan¬ 
tum perturbation theory. Plotting the energy levels as 
a function of some parameter, k for example, one can 
recognize that the mechanism of the enhancement due to 
RAT is similar to CAT: as the parameter is varied, the 
states forming the reference doublet, |T)/) in the present 
case, come close to a third state. They interact with each 
other, and in the interaction regime the splitting between 
the reference doublet becomes large, resulting in a spike 
fl3 |. Note, however, that the staircase structure formed 
with the plateau and steeply decaying regime has never 
been found at least in the completely integrable systems 
studied so far. 


III. STAIRCASE STRUCTURE WITH 
RESONANCE SPIKES 

A. Resonance spikes and the third states 

Each spike observed in Fig. 0b) appears as a result 
of energy resonance between the doublet |T)/) and a cer¬ 
tain third state. The spikes mostly appear in the plateau 
regime, but sometimes they are situated in the steeply 
decaying regime. In Fig. 0 we first demonstrate which 
type of third states are actually involved in the creation 
of spikes. In the original framework of the RAT theory, 
the predicted spikes are associated with the states sup¬ 
ported by nonlinear resonances in the inner torus region, 
encircling central elliptic fixed points, ( q,p) = (±1/2,0) 
in the present case. However, two of spikes in the first 
plateau appear as a result of resonance with the states 
associated with an outer transversal torus and the spike 
located at the end of the first plateau h = 1/27 is as¬ 
sociated with the state localized on the unstable fixed 
points (q,p) = (0,0) and (—1,0) (see Fig. 0b)). This is 
not surprising since the present eigenstates are Floquet 
states, so the eigenphase rE n /h of Eq. 0 can satisfy 
the resonance conditions E n — Eg = mh/r (n, £, m€Z). 
Therefore the quasi-energies of our reference doublet can 
resonate with a state associated with an outer transversal 
tours. Such situations are out of the scope of the theory 
of RAT, but as will closely be discussed in section IIV1 
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FIG. 3: (Color online) Husimi representation of the third 
states for r = 1, which resonantly interact with the refer¬ 
ence doublet in the first plateau regime: (a) 1/h = 18 (at the 
middle of the first plateau), (b) h = 1/27 (at the end of first 
plateau), (c) and (d) h = 1/44 (at the second plate regime). 
In (a), (c) and (d) we only show |4>J) states out of each tun¬ 
neling doublet. In (b), the state I'Fj) is not the one forming 
a doublet as mentioned in [23j|. Upper left box represents the 
size of effective Planck’s cell. 


the outer torus states play a crucial role in the formation 
mechanism of the staircase structure. 


Figures @|a) and (b) demonstrate the splitting A E 0 
(in the back panel), together with the behavior of dou¬ 
blet and the third state energies (in the floor panel) as 
a function of the parameter k. When a spike appears in 
the AEq vs 1/h plot, there always exist spikes in the plot 
of AEq vs k nearby. However, even if a spike appears in 
the AEq vs 1/h plot, it does not necessarily mean that 
one exactly hits a spike in the plot of AEq vs k. These 
figures reveal that the third state for 1/h = 27 (at the 
end of the first plateau) is the 10-th excited state 
and this state is, as shown in Fig. [3] (b), localized on the 
unstable fixed point (q,p) = (0,0). On the other hand, 
1/h = 44 (in the middle of the second plateau), there ap¬ 
pear two spikes in the range under observation, one is the 
doublet composed of the 8-th excited states whose 
symmetric state is shown in Fig. [3] (c), and the other 
is also given as a doublet of excited states, |4 , 2 q) whose 
symmetric state is shown in Fig. [3] (d). As demonstrated 
in Figs. [3] (c) and (d), both doublets |\l/g ) and |\I/± 0 ) are 
supported by elliptic inner and transversal outer KAM 
curves, respectively. 


B. Absorbing operator 


Since the present situation, the case of r = 1, is not 
far from the integrable limit, eigenstates \ Jg) in the in¬ 
tegrate limit well approximate eigenstates |\I/ n ), i.e., 
(Jel'&n) ~ 1- On the basis of this observation, we in¬ 
troduce the following absorbing operator [S3, [25| 

# = i -?Eiwi- (5) 

z leL 

Here T < 2 represents the absorbing strength, and 1 the 
identity operator. The summation runs over a given in¬ 
dex set L, which we choose appropriately depending on 
which states we want to suppress [34j . Below we con¬ 
sider the right eigenvalue problem for the absorbed (non¬ 
unitary) evolution operator 

U°\* n ) = u n \* n ), (6) 

where 

U° = PU, (7) 

and u n = e ~ lE n.T/h . f 0 ]l 0w i n g argument holds even 

if one considers the left eigenvalue problem. 

First we will discuss what we can expect in perturba¬ 
tion theory with respect to the absorbing strength T. It 
is easy to show that a standard perturbative calculation 
up to the second order provides 

U n — u n • z n , (8&) 


where 


z n 


— i _ 77 E i a b"i 2 + 

Z iGL 


r 2 

tEE 

l€/L m^n 


\ a l,n a e,m\ 2 
'Un 'Um 




(8b) 


and 


ri) • (8c) 

The right (absorbed) eigenstate is also given as 

|^n) - l^n) - W E Bm ’ n \®m) (9a) 

mzfcn 
r N -1 

= l^ro) ^ ^ ( 73. m n Q,] cr n\Jk)^ (9b) 

k =0 m^n 

where 

B m ,n = £ (9c) 

ttL Un - Um 

For 0 < T <2, (absorbed) quasi-energies E n are no 
more real because the second-order term in the pertur¬ 
bation expansion Eq. (18bl) becomes complex while the 
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(a) (b) 




FIG. 4: (Color online) Energy splittings AEo (black lines in the back panel) and the reference doublet (red lines in the floor 
panel) and the related third state energies (blue and light blue lines in the floor panel) as a function of the parameter k. For 
each curve, the corresponding energy level is put in the figures. Gray lines indicate energies of the states irrelevant to creating 
spikes. In the upper panels, no absorber was applied for (a) 1/h = 27 and (b) 1/h = 44. In the lower panels, the strength of 
the absorber is set as T = 1 for (c) 1/h = 27 and (d) 1/h = 44. The index set L is given as (a) L = {10} and (b) L = {8}. 


first-order term is real-valued. The eigenvalue u n is then 
shifted as 

arg u n = arg u n + arg z n . (10) 

By applying the absorber, the coupling between the 
absorbed and the rest of eigenstates is suppressed. This 
could be regarded as an inverse procedure of what is 
done in typical perturbation theory such as RAT the¬ 
ory, in which one starts with some unperturbed states 
| J n ) and build up desired eigenstates |T„) by adding per¬ 
turbation terms. The present absorbing method is, in a 
sense, to subtract perturbed terms | J n ) from the final 
state IH/n). Therefore, applying the absorber in this way 
would be a test to check whether the final state |\F n ) 
could be obtained as a result perturbation in terms of 
unperturbed states | J n ), and, if so, which unperturbed 
states are involved in the perturbation procedure. The 
present absorbing method is equivalent to the one used 
in the open quantum systems, e.g., (TTI. [26j - [2?l . in which 
the absorbers are adopted as the Heaviside step function 
(x\J) = H(x) or the Dirac delta function (x\J) = S(x). 

The efficiency of the absorbing method is demon¬ 
strated in Figs. He) and (d). For absorbed quasi¬ 
energies, the tunnel splitting A E± is defined as A E n = 
E~ — E //, however we note that E± has an imaginary 
part when T > 0. Each corresponds to the case where the 
absorber with F = 1 is applied to the case shown in Fig. 
He) and (d), respectively. Here the index set L is chosen 


as L = {10} in the case of 1/h = 27, and L = {8} whose 
member corresponds to the doublet of the symmetric and 
anti-symmetric state | Jg ) for l/h = 44. Here {£} repre¬ 
sents | Jp), and the member \jf) in the index set L is 
chosen in such a way that it maximally overlaps with the 
third state that is interacting the reference doublet IT^) 
and responsible for creating the spike. 

As clearly shown, energies of the associated third states 
gain certain amount of the imaginary part and pushed 
out to the complex plane, resulting in vanishing the 
spikes. The effect to the other states is almost negli¬ 
gible. However, as seen in Fig. Hd), the right-hand peak 
with shorter height still remains since we have not in¬ 
clude the states | J^) in the absorber. As mentioned in 
the end of the previous subsection IIII Ai there are two 
sets of doublets which are involved in avoided crossings 
in question. 


C. Staircase structure 

In the previous subsection, we have selected out ab¬ 
sorbing states by plotting energy levels around each 
avoided crossing and then judging by hand which states 
should be included in the set L , that is, it was necessary 
to refer to the figures like Fig. H We now introduce a 
systematic procedure to choose the absorbing states nec¬ 
essary to suppress the observed spikes. 
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FIG. 5: (Color online) Tunneling splittings AEo for the evo¬ 
lution operator U° are plotted as a function of 1 /h. The range 
of interaction in increased as (a) s = 1 , (b) s = 2, and (c) 
s = 3. Black curves show the unitary case (r = 0) and each 
colored curve the case with absorber, and the colors distin¬ 
guish the absorbing strength T (see the right-hand color bar). 
The splitting AEo in the integrable limit is shown as black 
dashed lines in each figure. 


shown in Fig. 0Jb), some spikes come close to each other 
in the second plateau and a single absorber is not enough 
to suppress the interaction with the reference doublet. In 
the case presented in Fig. QJb), the third state respon¬ 
sible for the left-hand peak is the state supported by an 
elliptic torus inside the KAM region and the right-hand 
one is supported by a transversal torus. As we further 
add the corresponding absorbers in this way, the peaks 
surviving in the s = 2 case gradually disappear, and the 
curve almost converges at s = 2. 

It would be worth emphasizing that steeply decaying 
regions are not affected and robust against the the ab¬ 
sorber applied on plateaus. This strongly suggests that 
the influence of spikes is well localized in each plateau, 
not like the situation suggested in El- This observation 
also supports our hypothesis; the splitting curve should 
be viewed as a staircase structure accompanied by spikes, 
not as spikes bringing the staircase. 


IV. MECHANISM GENERATING THE 
STAIRCASE STRUCTURE 


The most natural criterion to achieve this would be 
to check the energy difference from the reference dou¬ 
blet: d(E±) = | Eq — E ±| (n = 1,2 ,•••), because the 
spikes appear when the reference doublet and a certain 
third state are energetically close to each other and form 
avoided crossings. We rearrange the states in as¬ 
cending order of d(E±), and the corresponding integrable 
base \Jn) as well. The one-to-one correspondence be¬ 
tween |’F^ : ) and | J„ ) is again ensured since the condition 
|(J±|’F±}| ~ 1 is now satisfied for r = 1. Then the set 
of absorbing states containing the first s doublets in the 
sense of the energy distance reads 

L s = {1, 2, ■ ■ • , s}, (11) 

where we drop from the list of L s the states not belonging 
to the same parity as |Tq ). We must recall that the 
ground state |T^) has the symmetry with respect to the 
translation in addition to the mirror transformation. 

Note that the set L s of absorbing states depends on 
the value of h, so has to be determined for each h. As 
explained below, the reason why we consider the cases 
s > 1 is that the energetically nearest state L\ from the 
reference doublet is not sometimes sufficient for killing 
the coupling with the reference doublet. 

Figure [5] plots the splitting of quasi-energy Eq eval¬ 
uated for the operator ( 7 ° as a function of 1/h. In the 
s = 1 case, we see that some spikes, especially in the first 
plateau, disappear with increase in T. As shown in Fig. 
[2 the absorber pushes the third level into the complex 
domain, and the coupling with the reference doublet is 
suppressed. 

However, we notice that some spikes still remain in the 
second plateau regime. This is due to the fact that, as 


The main message in the previous section is that the 
staircase-shaped skeleton is formed in the splitting curves 
and spikes are superposed on it. In this sense we may 
say that the origin of the enhancement of the tunneling 
probability traces back to such a staircase structure. In 
this section, we study the mechanism creating the stair¬ 
case structure by introducing the renormalized basis, and 
show the reason why this only appears in nonintegrable 
systems. 


A. Instanton-noninstanton transition 

As shown in (2lj| . the TNI transition could be well 
captured by renormalized perturbation theory. An im¬ 
portant finding there was that a remarkable quenching 
of renormalized transition matrix elements explains the 
I-NI transition. In particular, without using highly renor¬ 
malized integrable Hamiltonian as unperturbative bases 
one could not identify the mechanism behind the transi¬ 
tion. 

For this reason, we also apply the same perturbation 
scheme to pursue the origin of the staircase structure. In 
essence, renormalized perturbation theory makes use of 
the Baker-Campbell-Hausdorff (BCH) expansion [2lU29l |: 

no'i 

where 

M 

(t p) = #1 (w P) + 

j=3 

(jGodd int.) 




eff ? 






U « Um = exp 


(13) 














1/h 


FIG. 6: (Color online) (a) Left panel shows the error 1 — 
|'Fn)| 2 as a function of the BCH order M in the case of 
h = 1/63. Not only the grand state, n = 0, but also excited 
states up to n = 74 are examined, (bl) The amplitude at 
q = 0 of the ground state obtained by the perturbation cal¬ 
culation, and (b2) the amplitude of the maximal mode of the 
contribution spectrum (see the text). The exact amplitude 
for the ground state |'I f o) at q = 0 is shown as blue curves in 
(bl) and (b2). The solid black dashed and dotted curves show 
the exact level splitting AEq and the level splitting AAq M ' ) 
of the integrable basis for reference, respectively. 


Here Hj denotes the j-th order term in the BCH series. 
Explicit forms for the first few terms are found as 

Hi(q,p) = T{q) + V(p), (14a) 

H 3 (q,p) = ^{[T, [T, V]] - [V, [V,T]]), (14b) 


where the terms Hj for even j are equal to zero thanks to 
the symmetrized form of U. The first order BCH Hamil¬ 
tonian is identical to the continuous time Hamil¬ 
tonian and higher order BCH Hamiltonians H^ 1 ' 1 are 
expressed as nested commutators. We denote the eigen¬ 
functions of the integrable Hamiltonian H by | 

hW)\jM) =E M\j[ m) ). (15) 

We first check the validity and efficiency of renor¬ 
malized perturbation bases by examining the error 1 — 
|(Jn M '|'k„)| 2 of the approximation. As shown in Fig. 
[HKa), the BCH states becomes better approximation to 
the corresponding eigenstate as the expansion or¬ 
der M increases. Note also that the expansion works for 
the lower energy eigenstates as compared to the higher 


excited states. This is, however, not a convergent expan¬ 
sion: the error 1 — |\E'„)| 2 starts to grow when the 

expansion order M exceeds a certain optimal order. 

Such highly efficient integrable approximation ensures 
the validity of renormalized perturbation, in which the 
difference A Um = U — Um could be regarded as a per¬ 
turbation (2lj . As also shown in Fig. (HKbl), the results of 
the 1st order perturbation calculation are in an excellent 
agreement with the exact ones, and even the staircase 
structure could be reproduced. However we would like 
to remark that although perturbation theory, not neces¬ 
sarily the present one, works well, this does not tell us 
anything about the underlying mechanism generating the 
staircase. 

As shown in Fig. [6fbl), the splitting A E is strongly 
correlated with the amplitude of the eigenstate at q = 0, 
and characteristic patterns appear around q = 0. As seen 
in Fig. [3(a), the eigenstate |^o) for r = 1 in the instanton 
regime is, as expected, well fitted by the one | Jq M ^) in the 
integrable bases, whereas the integrable approximation 
does not work any more and further structures appear 
in other regions. In the first and second plateau, the 
curve bends in a convex way (see Fig. (3b) and (d)), 
but in the first steeply decaying region the curve bends 
in a downward direction at q = 0 and takes a concave 
structure (see Fig. 13c)). 
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FIG. 7: (Color online) The black curve in each figure shows 
the eigenstate |'Fo) for r = 1 in the (a) instanton, (b) first 
plateau, (c) second decay and (d) second plateau regime, re¬ 
spectively. The dashed curve displays the integrable eigen¬ 
state (q\Jf) AI ^) at the corresponding h value, and colored ones 
the integrable components (q\Jg M ' > ){Jg M ' > \'i’o) at q = 0, where 
the value of £ is put in each figure. Note that the structure 
around q = 0 can be well reproduced by the maximal mode(s) 
of the contribution spectrum (see the text). 
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To explore the nature of wavefunctions at q = 0, we 
here introduce a spectrum decomposition at each position 
q in terms of integrable bases 

N-l 

<g|^o> = Coi 4 m) ( 9) ( i6 ) 

e=o 


where 


Con ( e M) (q) = (q\jj M) )(jj M) \V+). (17) 


We call such a decomposition the contribution spectrum 


21]. In the following discussion, we focus only on the 


symmetric ground state IH/g"). As we mentioned in Sec. 
mi each eigenstate has a symmetry with respect to the 
mirror transformation II g and the translation T r . There¬ 
fore, the basis | j\ M ' > ) that has the same symmetry as 


|\I/[j") is only used for the contribution spectrum. 

As shown in Fig. [Htb2), since the maximal mode of 
the contribution spectrum quite efficiently describes the 
behavior of the splitting, we can deduce that the staircase 
structure must be characterized by the maximal mode. 
Indeed, in Ref. [21], we have shown that the instanton- 
noninstanton (I-NI) transition could be explained as the 
switching behavior of the most dominant component in 
the contribution spectrum; from the one representing the 
instanton contribution to broad components supported 
around the separatrix of a central unstable fixed point. 
Below, we present that the dominant component controls 
not only the transition from instanton to noninstanton 
but overall signatures in the staircase structure. We will 
explain this by showing contribution spectra for several 
values of 1 /h, which are presented in Fig. [5] 

First of all, as mentioned just above, we notice that the 
contribution spectrum is mainly composed of two peaks 
with distinct characteristics. The first one is a sharp 
peak located at E = Eq M \ and the second is composed 
of many components, whose center is situated around the 
separatrix energy. A small peak sometimes appears on 
the broadly spread components as a result of the interac¬ 
tion with a third state. The first sharp peak at E = E q M ^ 
originates from the instantion contribution that has a 
maximal overlap with the ground state |'F]j"), so we here¬ 
after call it the instanton peak. We stress again that the 
instanton peak at E = Eq M ^ can be recognized only when 
we prepare higher order BCH expansions (M = 7 for the 
present calculation), otherwise the instanton peak is not 
isolated from the others and could not be identified. 

In the instanton decay (the first steeply decaying) 
regime, which is seen in the case of l/h = 12 in Fig. 
[Sfa), the instanton peak dominates the other compo¬ 
nents. As a result the amplitude of the ground state | d'[j") 
at q = 0 is well described by the integrable Hamiltonian 
base | Hence the instanton behavior should and 

is actually observed. 

With increase in l/h, the height of both peaks, the in¬ 
stanton peak and the broad components centered around 



FIG. 8: (Color online) The contribution spectrum Con^ 7 - 1 
(in logj q scale) at q = 0 is plotted as a function of the energy 
Ep M> . The BCH order M = 7 was used. The values of h 
are indicated in each figure. Each panel respectively shows 
the case (a) before (yellow) and after (red) I-NI transition, 
(b) before (red) and after (blue) the transition from the first 
plateau to the second decaying, (c) before (blue) and after 
(red) the transition from the second decaying to the second 
plateau, (d) before (red) and after (blue) the transition from 
the second plateau to third steep decaying regime. The dot 
represents the maximal mode in each spectrum. We have used 
a yellow-colored curve in the region where the maximal mode 
is given by the instantion contribution, a red-colored when 
the maximal mode energy is above the separatrix energy, and 
a blue-colored below the separatrix energy. The thick solid 
line represents the separatrix energy and red dotted lines the 
energies satisfying the condition E = Eq AI ' > + mh/r (m = 
0,1, 2, • • ■) (see the discussion in subsection IIV Bl) . 
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the separatrix, gradually drop, but the speed of the for¬ 
mer is much higher than that of the latter, eventually 
resulting in the switching of the role of the dominant 
contributor from the instanton to the top of broad com¬ 
ponents. An important remark is that the support of the 
state associated with the top of the broad components is 
outside the separatrix, meaning that the ground state is 
most dominantly coupled with an outside state [U]. We 
notice in Fig. [6fb) that, exactly at this switching mo¬ 
ment, the first instanton decay turns to the first plateau, 
and eigenstates show the convex structure around q = 0 
(see Fig. Ob)). In the perturbation calculation, it is 
also crucial to include outer torus states into unperturbed 
bases to reproduce the convex structure at q = 0, oth¬ 
erwise the resulting wavefunction cannot bend upward 
at q = 0 since it is merely a superposition of exponen¬ 
tially decaying states. It is important to note that not 
only such a convex structure just after the transition but 
also neighboring structures around q = 0 could be well 
reproduced only by the maximal mode in the spectrum 
Con^ M ^ (see Fig. [HIb2)). 

In any case, the maximal mode in the spectrum 
Conl M ^ can be a good indicator for the value of eigen¬ 
states at q = 0 and thereby the splitting A E. The max¬ 
imal mode in the contribution spectrum Con[ M \ which 
is shown using color-coded dots in the Fig. [9] well traces 
the staircase structure of the exact splitting A E n , and 
the value of eigenstates |(g|\E' n )|^ =0 at q = 0 as well. We 
will fully make use of this fact hereafter. 

As we further increase in l/h, the instanton peak is 
completely overtaken by the broadly spread components 
(see Fig. HKb)) and this ordering is fixed and never turned 
over. We also emphasize that the estimation of the crit¬ 
ical Planck’s constant h c at which the I-NI transition 
occurs becomes a bit imprecise if we use the lower order 
BCH series. 

As we increase l/h after the I-NI transition, the sup¬ 
port for the maximal mode of the contribution spectrum 
further approaches the separatrix, which is shown as in 
Fig. Htb), and eventually it goes into the inner tours re¬ 
gion in excess of the separatrix. At this moment, we re¬ 
alize that the splitting curve changes the behavior from 
the first plateau to the second steeply decaying regime 
(see Fig. [Htb)). At the same time, the structure of eigen¬ 
states at q = 0 changes from the convex to concave shape 
(see Fig. [He)). 

With further increase in l/h the maximal mode also 
shifts to the left. On the other hand, another peak is born 
at the right-hand edge of broad components, and now the 
competition comes into issue between the those peaks, 
the one playing a major role in the I-NI transition, and 
the new one at the right-hand edge. As noticed in Fig. 
Etc), the switching of the dominant contributor again 
takes place between these two peaks, and at this moment 
the splitting curve turns from the second steeply decaying 
to the the second plateau regime. 

After such a transition, the overtaken peak, the one 
playing a role in the I-NI transition, is gradually ab- 



FIG. 9: (Color online) (a) Maximal modes plotted as a 

function of the inverse Planck’s constant l/h. We have used 
yellow-colored dots in the region where the maximal mode is 
given by the instanton contribution, red-colored ones when 
the maximal mode energy is above the separatrix energy, and 
blue-colored below the separatrix energy, (b) Maximal modes 
for the absorbed eigenstates |4'„). The absorption procedure 
is the same as that introduced in section Em Here we use the 
integrable bases | j| M ^) as the absorber, and the absorption 
parameters are chosen as s — 3 and P = 0.4. The rule for color 
coding is the same as in (a). In both calculations, the BCH 
order M = 7 was used. In (a) the exact eigenfunction 
at q = 0, integrable basis ( q\Jn M "*), and energy splitting AI?o 
are shown as solid, dotted and broken curves, respectively. 
In (b) the solid curve represents the absorbed eigenfunction 
(ql’I'o') at q = 0, and dotted and broken ones are the same as 
in (a). 


sorbed into the spectrum envelope. However it leaves 
a clear trace in wavefunction: As shown in Fig. (Hd), 
the shoulder or bulge observed in the neighboring region 
around q = 0 is well reproduced by the component that 
has played a role in the I-NI transition. The convex struc¬ 
ture observed in the first plateau is pushed outward by 
the newly born component, and then it appears as shoul¬ 
ders. In other words, the history of the staircase struc¬ 
ture in the splitting plot is properly recoded in the tail 
of wavefunction, not necessarily at q = 0. 
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The staircase structure in the splitting plot could 
therefore be explained by the successive switching pro¬ 
cess of maximal modes, and passing through the separa- 
trix, that is, whether the support of the maximal mode is 
inside or outside the separatrix. Figure [9] illustrates that 
the staircase structure of the splitting curve can be un¬ 
derstood by the position of the maximal mode: whether 
its support is outside or inside the separatrix. 

We have verified, as shown in Fig. [9jb), that even if 
we suppress the peak standing on the broad peak com¬ 
ponents using the same absorber technique in subsection 
IIIIB1 this switching process still survive. This implies 
that the switching does not occur specifically between 
the resonance peaks appearing in the contribution spec¬ 
trum, but overall deformation of broad peak components 
controls it. We could identify at least the third and fourth 
transition and confirmed the same scenario applies. 


B. Anomaly of eigenfunctions in the action 
representation 

As shown above, we could attribute the emergence of 
the staircase structure to the successive switching of the 
dominant component in the contribution spectrum. In 
this subsection, we explain why the quantum number of 
the dominant component gradually shifts with increase 
in 1 / ft,, passing through the separatrix, and also explain 
why this causes the change in the slope of the splitting 
curve. 

For this purpose, we examine the behavior of the ex¬ 
pansion coefficient an d the integrable eigen¬ 

function (q\J^ M ' > ) at q = 0 separately. Note that the 
product of these two terms constitutes each component 
in the contribution spectrum Con^ M -*. We here call 
<4 M) |1/o> the eigenfunction in the action representation. 

First of all we remark that the value h/r becomes a 
fundamental energy unit in our system. This is because 
the present system is driven by the periodic force with 
period r, so t/t x h = h/r becomes a fundamental 
energy unit, and the energies specified as = E^ 1 ' 1 + 
mh/r (in = 0, 1 , 2 , • • •) may invoke quantum mechanical 
resonances. In Fig. [8] we have shown such energies as 
dotted red lines. In the following, we first describe a 
signature of {q\J^ AI ' > ) at q = 0 and then discuss anomaly 

found in (J» M IHfo}- Combining these, we finally explain 
the mechanism of successive switching in the contribution 
spectrum. 

As shown in Fig. [lUl the amplitude of the integrable 
eigenfunction (q\J^ M ^) at q = 0 shows exponential depen¬ 
dence on the energy E as far as the energy is less than 
that of the separatrix (left side of the thick black line in 
Fig. [TUI). When plotting ( q\j[ M ■*) with a fixed energy 
one also finds exponential decay as a function of 1/h (see 
Fig. Ob)). This is an expected behavior since {q\J^ M ^) 
is just an eigenfunction of an integrable Hamiltonian, no 



FIG. 10: (Color online) Amplitude of the integrable eigen¬ 
function (q\Jf Al ^) at q = 0. The 7th order BCH Hamiltonian 
was used. The black vertical line shows the energy of separa¬ 
trix. 


matter large the expansion order M is. 

On the other hand, above the separatrix energy (right 
side of the thick black line), we see that the amplitude of 
the integrable eigenfunction (g| keeps almost con¬ 

stant. This is also reasonable because each {q\j[ M ' > ) has 
its supports on a transversal invariant torus outside the 
separatrix, so the connection is not made via tunneling 
but real classical processes, thus resulting in no decay as 
a function of the energy. 

In contrast, the nature of the eigenfunction |\ko) 

in the action representation is highly nontrivial. As 
shown in Fig. fTlT a). there exists a sharp peak at Eq M \ 
which represents the instanton contribution, and then 
the value of suddenly drops to reach a small 

level. Then it forms a non-decaying region in which the 
value of {J[ Ml | T 0 ) does not decrease, rather increases 
gradually until a small peak which is close to the en¬ 
ergy which is specified by the relation E = E ^ 1 ' 1 + h/r 
[3u . f32l ] . This peak originates from the resonance of the 
associated states with the periodic forcing inherent in our 
model. The non-decaying region means that as long as 
the eigenphase difference is less than h/r the contribu¬ 
tion from the associated integrable basis states is almost 
equal. It is beyond this resonance that the exponen¬ 
tial decay common in the ordinary tunneling tail takes 
place. As presented in Figs. [TiTbWc). the presence of 
non-decaying region is not limited to the ground state 
but appears in exited states as well. Also note that over¬ 
all features are reproduced by just one-step time evolved 
wavefunciton which is expressed as (J^ M ' > \AUM\Jn M ' > )- 
The latter is consistent with the observation that pertur¬ 
bation theory based on the BCH basis works well (see 
Fig. [U])b)). We emphasize that these are all observed 
only when the order M of the BCH approximation is 
large enough and also universally appear in the eigen- 
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function of quantum maps [U H(| HU • 

A particularly important fact is, as shown in Fig 1121 
that the decay rate of the height of the non-decaying re¬ 
gion as a function of 1/h is extremely slow, as compared 
to the region E^ M ' > > h/r. This clearly distinguishes and 
characterizes the two regions, below and above the reso¬ 
nance energy E = E^ 1 ' 1 +h/r. We should make clear the 
underlying reason behind the observed power law decay 
in both regions, but the observed energy, so the corre¬ 
sponding classical structure as well, moves with increase 
in 1/h in the current setting, which makes difficult to 
apply a straightforward semiclassical argument. 

In addition to the resonance peak at E = E q M ^ + h/r, 
a sequence of peaks implying the higher order reso¬ 
nances appear at E = E^ I ' > + mh/r (to integer) (see 
Fig. Illl all. In conjunction with resonance peaks, there 
also exist narrow non-decaying regions just below each 
peak as the non-decaying region appearing in the region 
E—Eq AI ^ < h/r. Such a sequence of non-decaying region 
is not so sharply identified in Fig. EH a), but it becomes 
clearly visible as we increase 1/h. We can therefore divide 
each sector E^ 1 ' 1 +mh/r < E < E^ + (m+l)h/r into 
two characteristic regions; the one showing faster decay 
with 1/h and the other having quite slow decaying char¬ 
acter. A detailed explanation will be presented in our 
forthcoming paper [321] . and we just show in Fig. 1121 the 
difference of the decay rate by measuring it in the middle 
energy in each sector. As is seen, the decay rate in the 
region Eq^ + h/r < E < Eq^ + 2h/r is much slower 

than in the next sector E^ 1 ' 1 +2h/r < E < E'q M ' > +3 h/r. 
Although we do not specify in which characteristic region 
the middle point energies used to measure the decay rate 
is contained, it is enough, in the following argument, to 
notice that the decay rate much differs in each sector. 
Also note that such resonance peaks with the same na¬ 
ture also appear in exited states as also shown in Fig. 
EHb)-(c). 

Putting all the pieces together, we can now understand 
why successive switching in the contribution spectrum 
generates the staircase structure. In the first decaying 
(instanton) region, the instanton is the most dominant 
and broadly spread components provide only negligible 
contributions, as explained in the previous subsection. 
The height of the instanton peak decays exponentially 
with 1/h as expected. However, in this region, the largest 
component in the broadly spread components is outside 
the separatrix (see the yellow curve in Fig. fTUI aH. mean¬ 
ing that the separatrix energy is contained in the non 
decaying region of (|’P 0 ). Since (q\jj, M ' > ) keeps con¬ 
stant when the position q is outside the separatrix and 
the decaying speed of (jj M ^ |H/o) is so slow as shown in 

Fig. [12] its product Con^ M) also decays much slower 
than the instantion peak. Thus, at a certain critical l/h c , 
the instanton component is overtaken by the dominant 
component in the broadly spread components. This is 
nothing but the I-NI transition [gif . 



FIG. 11: (Color online) Eigenstates |'Fn) in the action 

representation plotted as a function of E for h = 1/80 for 
(a) n = 10, (b) n = 15 and (c) n = 25, respectively. The black 
curves show the matrix elements (j/ M ' l \AU\j/ M ' > ). Here we 
used the 7-th order BCH Hamiltonian as the basis \J < 'p M> ). 
The black solid line and dotted lines respectively show the 
separatrix energy, and the energies satisfying the condition 

E = E ( 0 M) +mh/r (m = 0,1,2, •••)■ 


After the I-NI transition, as long as the position of 
the dominant contribution in the broadly spread com¬ 
ponents is outside the separatrix, the decaying behav¬ 
ior in the plateau of (T^ M ' , |'F 0 ) controls the product 

(g| Jg M) ITo). This explains the presence of plateau 

in the splitting curve. 
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FIG. 12: (Color online) The inverse Planck’s constant 1/h 
dependence of (jj M ^|'l/o). Difference of colors distinguishes 
the energy at which the value of |'Fo) is evaluated (see 
the inner panel). 

However, note that the position of the dominant com¬ 
ponent is determined by the edge of the plateau of 
(jj M ^ |^o), an d this edge is located around the value h/r. 
As a result, at a certain value of 1/h, the position of 
the dominant contribution passes through the separatrix 
(see Fig. [Sfb)). If such an event occurs, the separatrix 
energy is then situated in the region where (jj M) I'Fo) 
shows faster decay. This is exactly the moment when the 
splitting curve turns from the first plateau to the second 
steeply decaying region. 

The mechanism generating the next plateau is under¬ 
stood by observing ( j\ M ^ |’Fo) hr a wider range. As shown 
in Fig. HHa), a sequence of peaks appears at integer 
multiples of the fundamental energy unit h/r, and the 
decay rate of |\Ho) just below each resonance peak 

is again very slow as compared in the next sector, as 
demonstrated in Fig. [12] Hence the same switching pro¬ 
cess takes place repeatedly. We have actually checked 
that the mechanism explained here works at least until 
the third plateau, but we expect that this continues in 
larger 1/h regimes. 

In this way, we could explain the emergence of the 
staircase structure based on the nature of the action rep¬ 
resentation, which seems to be closely connected with the 
fundamental energy sequence whose unit is given as h/r. 
As was checked above, the fundamental energy sequence 
can induce quantum resonances, resulting in the spikes 
in the splitting curve. However, it should be noted that 
the appearance of quantum resonances is not a neces¬ 
sary condition for the presence of the staircase structure, 
as discussed in the subsection IIII Cl and IIV Al In other 
words, even if the resonance condition is not satisfied, 
a broadly spread or mild peak, whose width is almost 
comparable to the fundamental energy unit h/r, survives 
around the fundamental energy sequence. This is quite 


an anomalous situation because such a broad peak im¬ 
plies the existence of periodic oscillation of period r ac¬ 
companied by a rapid decaying process whose life time is 
comparable to the oscillation period itself [30| . 

We also characterize this anomaly from the view¬ 
point of semiclassical theory. If the leading-order 
semiclassical approximation works, the matrix ele¬ 
ment (j] M) |A17m| Jq M) ) should take a form of HI ~ 
^2 A 1 e~ lS ~</ h , where A 7 and S 7 respectively stand for 
the amplitude and classical action, and the sum 7 is 
taken over complex classical orbits satisfying given ini¬ 
tial and final conditions. In the semiclassical regime, we 
may neglect the h dependence in the amplitude A 7 , so 
the matrix element ( j[ M ^ \ AUm\Jq AI ' > ) is approximately 
expressed using the minimum imaginary action Im S~ /0 as 

~ g-ims 70 /?i. gj nce I m S- /0 is a purely classical quan¬ 
tity, the form h\n( \AUm\Jq M ' > ) should not depend 
on h. As will be shown in Fig. [T 6 ] this is indeed the case 
in the integrable system. On the other hand, Fig. [T3l 
shows that the matrix element \AUm\Jq M ^) does 

not follow the semiclassical ansatz in the non-decaying 
region, whereas the leading-order semiclassical predic¬ 
tion seems to work well beyond the non-decaying re¬ 
gion. Although it is necessary to check whether or 
not the leading-order semiclassical approximation indeed 
breaks in the non-decaying region, the observed sharp 
distinction would be an important signature characteriz¬ 
ing anomaly. 

According to these speculations, we are currently tak¬ 
ing two approaches to understand what was observed in 
the eigenstate (J^ A/ ^|H , o) in the action representation; 
one is a real semiclassical analysis which is based on 
the so-called classical-quantum correspondence principle. 
This could extract anomalous components hidden in clas¬ 
sical dynamics generated by the BCH Hamiltonian, and 
actually reproduce anomalous decay tails [30|. Another 
approach is to take into account higher-order effects in 
the semiclassical analysis. Since similar non-decaying or 
anomalous behaviors have been found in the model with 
discontinuity in phase space, observed phenomena might 
be liked to or have at least close similarity with diffrac¬ 
tion [3lj. This naturally leads us to the semiclassical 
treatment beyond the leading order. In any case, these 
are out of the scope of the present paper, and will be 
reported closely in our forthcoming paper |32|. 


V. SPLITTING CURVES IN INTEGRABLE 
SYSTEMS 

In the previous subsection, we discussed the underly¬ 
ing mechanism controlling the staircase structure of the 
splitting curve and found that anomalous tails in eigen¬ 
functions in the action representation play a key role. 
If such a feature is shared only in nonintegrable maps, 
we would not expect the enhancement of the tunneling 
probability in the completely integrable system. Below 
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FIG. 13: (Color online) Scaled wavefunction 

ftln \(Jg M ^\AU\Jf M ^)\ 2 as a function of Ep M ' > for several 
effective Planck’s constant h. The black solid line and 
dotted lines respectively show the separatrix energy, and the 
energies satisfying the condition E = Eq M ■* + h/r 

we shall explain, the nature of the splitting curve in the 
integrable system is totally different, although a seem¬ 
ingly common behavior is observed. 

For this purpose, let us consider the following classi¬ 
cally integrable Hamiltonian 

H(q,p) = H 0 (q,p)+£H 1 (q,p) (18) 

with 

H 0 (q,p) = (cos 2 q + cos 2 p)/2 + a(cos 2 q + cos 2 p) 2 , 

(19a) 

Hi (q,p) = cos 4 p — 6 cos 2 p cos 2 q + cos 4 q. (19b) 

This system was analyzed in 0 to examine the validity 
of RAT theory in a completely integrable situation. The 
authors have introduced a parameter (f> which controls 
the relative orientation of the classical resonance chains 
0 . Since the formulation of RAT theory do not take 
into account such orientation, RAT calculation could not 
follow the difference originating from it 0. 

As is seen from the Figs. fCT al) and (a2), the equi- 
energy surface has a local maximum between an unsta¬ 
ble fixed point ( q,p) = (0,0) and a stable fixed point 
( q,p) = (±^,0). Some equi-energy surfaces in the in¬ 
ner well, which appear around the stable fixed points 
(q,p) = (±f,0), have the same energies as those in the 
outer region. For e > 0, a classical nonlinear resonance 
chain is developed along the ridge between the inner well 
and outer region. 

We impose the periodic boundary condition on the re¬ 
gion (q,p) € (—7T, 7r] x (0, 7r] , and solve the eigenvalue 
problem 

mmn)=E^t)- ( 20 ) 

We then consider the splitting A E 0 = Eq — Eq of the 
ground and first exited states, both localizing in the inner 



FIG. 14: (Color online) Phase space portrait for the Hamilto¬ 
nian (1181) for a = —0.55 and (al) e = 0 and (a2) e = 5/1000. 
The black curves show the energy contour whose energy value 
is close to the maximum one. The black box put in the upper 
right corner represents the size of effective Planck’s constant 
for h = 1/5. (b) The splitting AEo (in log 10 scale) as a 

function of 1/h in the cases of £ = 0 (black dashed line) and 
£ = 5/1000 (black solid line). Yellow ad blue dots represent 
the maximal mode of the contribution spectrum Con( M ^ at 
q = 0 for £ = 0 and for 0 < t/N <1/2, respectively. The gray 
line shows the slope of the splitting curve for £ = 5/1000. The 
inset is magnification of a small 1/h regime. 

well. Here we take the innermost state in the inner well 
as the ground state and arrange the eigenstates in the 
same order as the standard map. 

Figure [CT b) gives the splitting AEq as a function of 
1/h. For e > 0, the splitting decays exponentially ac¬ 
companied with periodic spikes. All the features have 
clearly been accounted for if one applies the semiclassi- 
cal method using complex paths [l9(. The spikes appear 
as a result of the energetic resonance between the states 
localized in the inner well and outer region. The coupling 
strength could be evaluated using the imaginary action 
of complex trajectories which bridge classical disjointed 
equi-energy surfaces. 

It would be worth mentioning that for e = 0 the con¬ 
dition H (g, p) = 0 can be factorized into 

cos 2 q + cos 2 p = 0, (21) 

and 

cos 2 q + cos 2 p = — l/2a. (22) 

This shows that the invariant curves specified by m 
and (1221) are not connected even in the complex plane, 
thus no tunneling connection between the inner and outer 
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l/N l/N 


FIG. 15: (Color online) (a) The amplitude of (<7|Tf) at q = 0 
(dashed curves) and the action representation (J n |®o) (solid 
curves) as a function of the normalized quantum number l/N. 
Right panels give the contribution spectrum Con[ M ^ (log 10 ) 
as a function of the normalized quantum number l/N for (bl) 
a small 1/h regime, and (b2) a semiclassical regime. 


regions exists even though both are the surfaces with the 
same energy [33]. As a result, the splitting AEq exhibits 
single exponential decay without spikes. 

On the other hand, with careful observation of the 
splitting curve for e > 0 (see Fig. Mb)), we notice that 
there exists a crossover from one slope to another. In a 
small 1/h regime, the slope can be well fitted by the one 
for £ = 0, whereas the best fit curve, colored in gray in 
Fig. 0b), shows another slope in the large 1/h regime. 

Such a crossover or the change of the slope of the 
splitting curve reminds us of the plateau discussed in 
the nonintegrable situation. However, the origin and the 
underlying mechanism entirely differs from the previous 
one. This can be confirmed again by examining the con¬ 
tribution spectrum. Here we use the eigenstates |J„) 
as the basis states for the contribution spectrum, where 
H 0 (q,p)\J n ) = Ei 0) \J n ). 

As shown in Fig. I15l bl). a switching from the instan- 
ton to another mode also occurs, like the standard map 
case. However, in the integrable case, the position of the 
peak sits at the same value of l/N and does not move 
even if the value of 1/h is changed, while remember that 
it depends on 1/h and shift leftwards in case of the stan¬ 
dard map. Note here that l/N can be identified with the 
action coordinate. The reason for the peak position being 
fixed is simple; the peak appears as a result of the cou¬ 
pling between inner and outer surface, which is expected 
to occur in the RAT scenario. Alternatively stated, the 
origin of coupling is purely classical. Furthermore, as 
shown in Fig. [T6l the leading-order semiclassical anstaz, 
which was discussed in the previous section, works quite 
well for the eigenfunction in the action representation. 
These results make a sharp contrast to the standard map 



FIG. 16: (Color online) Scaled wavefunciton filn |(J^|'I'J)| 2 as 
a function l/N for several effective Planck’s constant h, each 
of which is taken at the values of off resonance positions in 
FigdHJb). 

case. We can see in Fig. [T7] that the maximal mode in 
the contribution spectrum well reproduces the structure 
of eigenfunction around q = 0, and its support is exactly 
an invariant curve with the same energy as that of the 
ground state. 

The presence of the crossover admits a simple semi¬ 
classical interpretation. As discussed in [l9j], there exist 
two different complex paths with different imaginary ac¬ 
tions. One corresponds to the ordinary instanton path, 
which runs from one well to another directly and the 
other is the path bypassing the classical resonance chain. 
In the semiclassical regime, since the latter one has a 
smaller imaginary action. On the other hand, in a small 
£ regime, it can happen that the instanton contribution 
is larger than that from the bypassing one, in spite of the 
magnitude relation of imaginary actions. This is because 
the prefactor, more precisely the coupling amplitude due 
to tunneling, comes into play in a relatively small 1/h 
regime. The observed crossover would be understood by 
taking into account not only the imaginary action but 
the coupling amplitude. This argument suggests that, 
in a larger e regime, the coupling with bypassing path 
gets larger, and the crossover point disappears when the 
value of £ exceeds a certain threshold. Note, however 
that the splitting curve cannot form the staircase struc¬ 
ture since we have at most two possible complex paths, 
and the underlying mechanism generating spikes has a 
purely classical origin as stated above. 


VI. SUMMARY AND OUTLOOK 

The focus of the present paper was put on clarifying the 
origin of the enhancement of tunneling probability in the 
nearly integrable system. We here measured the tunnel¬ 
ing probability by observing tunneling splittings plotted 
as a function of the inverse Planck’s constant. Typical 
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FIG. 17: (Color online) (a) The eigenfunction I'Fo) in the q- 
representation for e = 5/1000 (blue), the eigenfunction |Jo) 
(dotted) and the maximal mode in the contribution spectrum 
at q = 0 (green), respectively, (b) The maximal mode state 
in the Husimi representation for h = 1/5. The black thick 
curves represent inner and outer invariant curves with the 
same energy E = Co¬ 


features of the splitting curve commonly observed in non- 
integrable quantum maps are the existence of spikes and 
persistent departure of the splitting curve from the one 
predicted by the instanton. 

So far these have been discussed in the framework of 
RAT theory, but here we took a different perspective: the 
splitting curve is composed of the staircase-shaped back¬ 
bone accompanied by spikes. We have observed that, by 
introducing the absorber composed of integrable bases, 
spikes could be selectively suppressed if states interacting 
with the reference doublet are absorbed. More precisely, 
we have shown that the interacting or third state could 
be decoupled from the reference doublet when we take an 
integrable state that has maximally overlap with the cor¬ 
responding eigenstate as the absorber. The eigenenergy 
of such eigenstates responsible for generating spikes are 
pushed out to the complex plane, and spikes disappear. 
Our observation was that even though all the third states 
which resonate energetically with the reference doublet 
are suppressed in such a way the staircase structure sur¬ 
vives. Note that the efficiency of the present absorber 
comes from the fact that the regime we consider is close 
enough to the integrable limit, otherwise the absorber 
may affect irrelevant states in an uncontrollable way. 

The result strongly suggests the existence of non¬ 
trivial broad interaction between the reference doublet 
and other states. This was indeed confirmed by intro¬ 
ducing renormalized Hamiltonian, which is constructed 
using the BCH expansion, and used as basis states by 
which the reference state is expanded. In particular, we 
focus on the contribution spectrum at the origin q = 0 
since the amplitude of eigenfunctions at the origin follows 
the behavior of tunneling splitting. Here the contribu¬ 
tion spectrum introduced in (2lj represents components 
of renormalized states in the reference state. 


The contribution spectrum analysis clearly revealed 
that, in addition to the self component representing the 
instanton, there certainly exists broad interaction, and 
the behavior of such broadly spread components controls 
the staircase structure in the splitting curve. There are 
two key ingredients to explain the emergence of the stair¬ 
case: one is the behavior of the most dominant state in 
broad components, the other is anomalous tail observed 
in the eigenfunction in the action representation. Note 
that the renormalized bases are crucially important to 
capture these features, otherwise one could not explain 
the existence of the staircase structure and the anoma¬ 
lous tail part in the action representation as well. 

The dominant contributor in the broadly spread com¬ 
ponents switches from one to another, which was ob¬ 
served in the contribution spectrum. Such a switching 
phenomenon is driven by and liked to the existence of the 
fundamental energy sequence, which is further enhanced 
when the quantum resonance between unperturbed sys¬ 
tem and the periodic driving occurs. 

The origin of anomalous part in the action representa¬ 
tion should be explored more closely, which will become 
a primary subject of our forthcoming papers. The semi- 
classical analysis based on the correspondence principle, 
in which not complex but real classical orbits are used 
as input information. This efficiently works and turns 
out to extract anomalous components in classical dy¬ 
namics of the BCH Hamiltonian [3Cj • The analogy with 
the system modeling the diffraction, together with some 
speculations on anomalous behaviors of caustics appear¬ 
ing in the semiclassical analysis will be another approach 
[32j ■ The latter suggests that observed phenomena in the 
eigenfunction in the action representation are beyond the 
leading semiclassical description. 

These two key characteristics are, by their very nature, 
absent in the completely integrable system. Therefore, 
one could predict that the staircase structure does not 
appear in the completely integrable. We have confirmed 
this for a normal form Hamiltonian system, for which the 
validity of RAT theory was recently investigated. We 
have shown that a sharp contrast exists between inte¬ 
grable and nonintegrable systems and verified that the 
dominant contributor in the contribution spectrum for 
the integrable system sits at the same position and does 
not move as in the nonintegrable case. The absence of 
the staircase structure could simply be interpreted by the 
fact that there exists a unique dominant complex path in 
the semiclassical regime. 

Finally we would like to emphasize the importance of 
observing wavefunctions in the whole range, not focusing 
only on the amplitude at a specific point, like the origin 
q = 0 in the present case. As discussed in subsection 
IIV A[ with increase in 1/h, the convex structure around 
the origin, appearing in the first plateau, is pushed out¬ 
ward and forms shoulders in both sides. The same pro¬ 
cess happens repeatedly as one further increases 1 /h, that 
is, similar shoulders appear one after another. In this 
sense, we can find the trace of the staircase of the split- 
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ting curve in the tail pattern of wavefunction. This is 
also true for wavefunction in the action representation. 
There exists a significant difference between inner and 
outer tunneling tail, and this exactly results in different 
slopes of the splitting curve and thus staircase skeleton. 
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